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Abstract. The bright eclipsing and bursting low-mass 
X-ray binary EXO 0748—676 has been observed at sev- 
eral occasions by XMM-Newton during the initial cali- 
bration and performance verification (CAL/PV) phase. 
We present here the results obtained from observations 
with the EPIC cameras. Apart from several type-I X-ray 
bursts, the source shows a high degree of variability with 
the presence of soft flares. The wide energy coverage and 
high sensitivity of XMM-Newton allows for the first time 
a detailed description of the spectral variability. 
The source is found to be the superposition of a central 
(~ 2 x10 s cm) Comptonized emission, most probably a 
corona surrounding the inner edge of an accretion disk, 
associated with a more extended (~ 3xl0 10 cm) thermal 
halo at a typical temperature of ~ 0.6 keV with an indi- 
cation of non-solar abundances. Most of the variations of 
the source can be accounted for by a variable absorption 
affecting only the central comptonized component and 
reaching up to Nr ~ 1.3xl0 23 cm -2 . The characteristics 
of the surrounding halo are found compatible with an ir- 
radiated atmosphere of an accretion disc which intercepts 
the central emission due to the system high inclination. 

Key words: Stars: individual (EXO 0748-676/UY Vol) 
- binaries: eclipsing - X-rays: stars - accretion: accretion 
disc 



1. Introduction 

The X-ray binary EXO 0748—676 is unique in showing all 
types of variability commonly seen in different low-mass 
X-ray binaries. The source was discovered in outburst by 
the EXOSAT satellite. The observations revealed the exis- 
tence of sharp eclipses, intensity dips as well as repetitive 
type I bursts which were used to infer the presence of a 
neutron star as the compact object and to estimate the 
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distance to ~ 10 kpc (Parmar et al. 1986 hereafter PA86, 
Gottwald et al. 1986). More recent observations with the 
Rossi X-ray Timing Explorer (RXTE) have also shown 
that the source exhibits variable (~ Hz and ~ kHz) quasi- 
periodic oscillations (Homan et al. 1999, Homan & van 
der Klis 2000). The length of the orbital period (3.82 hr), 
combined with the duration of the eclipses, allows the ge- 
ometry of the system to be constrained with a high ~ 75° 
inclination and a ~ 0.5 Mq companion. The presence of 
residual flux during eclipses and of dips in the orbital light 
curve makes EXO 0748—676 typical of the ADC (accretion 
disc corona) sources in which the central compact source 
is surrounded by an asymmetric disc with the presence 
of a bulge in the outer edge and an extended central hot 
region (White et al. 1982). This complex geometry has 
prevented so far the different components to be isolated 
unambiguously and contradictory interpretations of the 
spectral variability have been put forward. From EXOSAT 
data, Parmar et al. (1986) concluded to the presence of a 
cut-off power-law emission whose part was strongly ab- 
sorbed. However from ASCA observations, Church et al. 
(1998) described the spectrum as the superposition of a 
point source (~ 2 keV) blackbody from the neutron star 
and a partially covered Comptonized emission, originating 
from the ADC. 

The very rich variability of EXO 0748—676 makes it an 
ideal target for XMM-Newton. The high throughput and 
large energy coverage of the EPIC cameras provide for 
the first time a detailed picture of the system. We present 
here the results obtained with the EPIC MOS and PN 
cameras. More results on the RGS Grating observations 
of the source can be found in Cottam et al. (2001). 

2. Observations and data reduction 

EXO 0748-676 was observed by XMM-Newton at four oc- 
casions during the CAL/PV phase, in revolutions 50, 55, 
59 and 67 of the satellite. Details on the XMM-Newton 
mission and on the EPIC MOS and PN cameras can be 
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Table 1. Log of EPIC observations 



Instrument 


Mode 1 /Filter 


SR 2 


Start time 


T 

J- cxp 








(2000, UT) 


(ks) 


MOS1+2 


FF/Thin 


59 


Apr. 04, 17:18 


20.0 


MOS1+2 


FF/Medium 


67 


Apr. 21, 04:11 


16.9 


PN 


FF/Medium 


67 


Apr. 21, 03:59 


18.1 


PN 


LW/Medium 


67 


Apr. 21, 09:50 


3.2 


PN 


SW/Medium 


67 


Apr. 21, 15:22 


18.1 



(1) FF Full Frame, LW Large and SW small window mode 

(2) Satellite revolution 



found in Jansen et al. (2001), Turner ct al. (2001) and 
Striider et al. (2001). 

Scientific data from EPIC MOS were secured in revo- 
lutions 59 and 67 (see Table 1). The data, obtained with 
a 2.6 s time resolution, were processed with the standard 
analysis of the XMM Science Analysis System. This in- 
volves the removal of bad or hot pixels, and of the elec- 
tronic noise. For the MOS, X-ray events with pattern 
to 12 (similar to grades 0-4 in ASCA) were selected. The 
gain was adjusted with reference to the closest observa- 
tion in the Closed Calibration position. The source events 
were selected by applying a mask of radius 1.5' centered 
on the source which corresponds to ~ 98% of the energy 
of the typical point spread function of the telescope (Gon- 
doin et al. 2000). The source is bright and some effects of 
the pile-up of events might be expected in the FF mode. 
The average counting rate in each MOS camera, except for 
the bursts, was found to be between (7-10) event /reading 
frame which is at the threshold for significant pile-up (see 
Fig. 6 in Ballet 1999). The pile- up rate is expected to be 
less than ~ 0.8% so that a moderate effect should affect 
the most central part of the image. 

EPIC PN data are currently available from revolution 
67 of XMM-Newton in three different modes (Table 1). For 
light curves we selected the valid event-patterns (single, 
double, triple and quadruple; pattern 0-12) while for the 
spectral analysis only single-events (pattern 0) were used 
for which an advanced spectral response matrix exists. To 
avoid pile-up effects only the data from the fastest read- 
out mode (the SW) was used for the spectral analysis. 

3. The varying light curve of EXO 0748-676 

3.1. The EPIC MOS observations 

Light curves were produced by selecting the source events 
in two different soft S (0.5 - 2 keV) and hard H (2 - 10 
keV) energy ranges and adding the counts from the two 
MOS cameras. The background rate is found negligible 
with mean values of ~ 0.04 and ~ 0.08 counts s _1 , re- 
spectively in the soft and hard bands, illustrating the very 



clean images obtained by XMM-Newton. Figures 1 and 2 
show the extreme variability of the source for the two ob- 
servations. In each observation, an eclipse is clearly seen 
in the hard band as well as several type I bursts with 
counting rates reaching up to 25 counts s _1 . 

The characteristics of the two eclipses are similar to 
those reported from EXOSAT (PA86) and RXTE (Hertz 
et al. 1997). The present lack of exact absolute timing 
prevents a precise comparison with the source ephemeris. 
However, the use of the provisional start times listed in Ta- 
ble 1 gives observed mid eclipses times which fall within 
17 s (April 4) and 38 s (April 21) of the predicted times 
using the linear ephemeris of Hertz et al. (1997). More 
tight constraints on the period changes will be provided 
when precise absolute timings will be available. Due to the 
very low background rates, the residual counts in eclipse 
can be measured for the first time very accurately. In 
the hard band, the residual counts are (0.50±0.20) counts 
s _1 and (0.37±0.18) counts s _1 on April 4 and 21 re- 
spectively, which corresponds to (10±4)% and (7±3)% of 
the corresponding mean flux outside eclipse. Significantly, 
this ratio, in the time interval corresponding to the hard 
eclipse, is (87±25)% and (32±10)% in the soft band. The 
behaviour of the source is different in the two observations 
as seen in Fig. 1 and 2. 

On 4 April 2000, the variability is mainly confined 
to the hard band with the presence of two typical dips, 
around phase ~ 0.7 of the 3.82 hr cycle (where phase 
is defined as the mid-eclipse time). There is no evidence 
in the soft band of the eclipse nor of the dips. The hard- 
ness ratio, computed as (H-S)/(H+S), is also shown at 
the top of Fig. 1 and 2. A significant softening is seen dur- 
ing eclipse, bursts and dips. On 21 April 2000, the source 
appears in a quite different and more active state. An in- 
tense flaring activity is seen in the soft band with two 
major increases of flux arising around phase 0.2 and 0.6. 
The hardness ratio is marked by a significant softening 
during these flares. 

A similar activity around the same phases has been 
noted from previous ASCA observations though above 1 
keV and with somewhat different energy bands (Church 
et al. 1998). The lack of information at lower energy has 
made the interpretation uneasy in this case. The very wide 
spectral coverage of EPIC allows here a more straightfor- 
ward separation between the lowest energies (below 2 keV) 
where the absorption effects are dominant and the higher 
ones (above 2 keV) where they are negligible (see below). 
We note also that our definition of the hardness ratio is 
quite different to that adopted in the 1-10 keV range for 
the EXOSAT observations (PA86). 

3.2. The EPIC PN observations 

The EPIC PN observations of EXO 0748-676 in FF and 
SW modes on April 21 (see Table 1) were performed ac- 
cidentally at almost identical binary phases. The light 
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Fig. 1. The MOS hard H (2 - 
10 keV) (bottom) and soft S (0.5 

- 2 keV) (middle) light curves 
of EXO 0748-676 as observed by 
XMM-Newton on 2000 April 4. 
Time in hours from 2000 April 4, 0:0 
UT. The curves are source counts 
obtained from the two MOS cameras 
and binned to 10 s. In each case, the 
underlying curve is the background 
rate accumulated in a region of the 
image free of sources. The hardness 
ratio computed as (H-S)/H+S) is 
shown on the top, where the curve 
have been smoothed by a 5-points 
average. Note the near total eclipse 
at ~ 19:40 UT and the two dips at 

- 18:30 UT and - 22:20 UT. Burst 
peak counts are cut in the figure 




Time (hour) 



Fig. 2. The MOS light curves for 
2000 April 21, with the same scale 
as Fig. 1. Time in hours from 2000 
April 21, 0:0 UT. Note the two sig- 
nificant soft flares around 06:30 UT 
and 07:40 UT 



curves are qualitatively identical showing marked features 
like flaring activity and dips at similar phases. Since dur- 
ing the FF observation also the MOS cameras observed 
EXO 0748-676 (Fig. 2) we show only the EPIC PN light 
curves of the SW observation in Fig. 3. 

4. Spectral variability 

Source spectra were accumulated from the EPIC PN data, 
using a circular extraction region of radius 40" around the 
peak surface brightness. Background spectra were created 
in adjacent regions of the same size. In view of the large 
variability of the source four representative spectra were 
chosen by selecting events in time intervals correspond- 



ing to the first and second soft flaring episodes (Fl and 
F2), the between- flare or quiescent interval (Q) and the 
dip part (D) of the EPIC PN SW observation on April 21 
(see Fig. 3). The count rate during this observation was 
well below the pile- up limit of the SW mode of ~ 100 
counts s _1 . Spectra were binned to obtain at least 20 
counts per bin. The most recent response matrices (ver- 
sion Sept. 2000) were obtained from the hardware groups 
and the spectra analyzed using the XSPEC package. 

Descriptions with a power-law and either a blackbody 
or a multi-temperature blackbody accretion disc model, 
as suggested from previous observations, were found inad- 
equate (see Table 2 where reduced % 2 values are given for 
various spectral models). In particular none of the models 
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Fig. 3. EPIC PN light curves of 
EXO 0748-676 from the small win- 
dow mode observation on 2000 April 
21 in the 0.5 - 2.0 keV (top) and 2.0 
- 10.0 keV energy bands. Data bin- 
ning is 30 s and time in hours from 
2000 April 21, 0:0 UT. The back- 
ground contribution is between 0.05 
and 0.2 counts s _1 in the soft band 
and 0.05 - 0.5 counts s _1 in the hard 
band. Note the similarity of the light 
curves during different binary orbits 
(see Fig. 2), in particular the stable 
phasing of the flaring events 



represents the spectrum around 0.7 keV where a steep de- 
cline due to an emission line feature can be seen. This line 
can be attributed to the O-VIII ion. To model the soft part 
of the spectra below 2 keV where more lines are visible we 
therefore tried thermal emission models, keeping in mind 
that the models available in XSPEC were developed for 
a thin plasma in collisional equilibrium, while the plasma 
which likely exists in EXO 0748—676 is photo-ionized by 
hard X-ray emission from a central source (see below). 

Best fits in the range 0.15 - 10 keV were obtained with 
a combination of a power-law (PL) and a Raymond-Smith 
thermal component (RS), each attenuated by individual 
photo-electric absorption (Morrison & McCammon 1983). 
Elemental abundances were treated relative to solar (An- 
ders & Grevesse 1989). It was checked that the fits with 
the Mewe et al (1985) instead of the RS thermal model 
gives fully consistent results including for the abundances. 

The four representative spectra were jointly fitted by 
the same model. In the power-law component, intensity 
and index were each kept as single free parameter tied to- 
gether for all four spectra. Only the column density atten- 
uating this component was allowed to vary individually. 
Similarly in the thermal component, the temperature and 
the elemental abundances were tied, but intensity and ab- 
sorption allowed to vary separately. The abundances of 
Si and higher Z elements were fixed to zero as is sug- 
gested by the RGS spectrum where no emission lines of 
these elements are found (Cottam et al. 2001). The other 
abundances were allowed to vary freely. The best fit RS 
model results in non-solar abundances of N (10±2), O 
(0.6±0.1), Ne (0.02±0.02) and Mg (0.31±0.23). The Mg 
line was only seen in the spectrum between flares. Table 
3 summarizes the best fit continuum parameters for the 
different intensity levels of the source. All quoted errors 



Table 2. Model fits to the flare spectrum (Fl) 



Model 1 


Xr/d0i 


powerlaw*wabs 
(powerlaw+bbody)*wabs 
powerlaw*wabs+bbody*wabs 
powerlaw*wabs+diskbb*wabs 


6.66/1071 
2.10/1069 
1.57/1068 
1.62/1068 



(1) Model components as defined in XSPEC 



were determined at A^ 2 of +2.7 around the \ 2 minimum. 
The spectra are shown in Fig. 4 with superimposed best fit 
spectral models. Residual features around emission lines 
seen below 1 keV are caused by calibration problems of the 
charge transfer losses in SW mode. In the current version 
of charge transfer loss correction in SW the event energies 
are over-corrected by about 10 - 20 eV at energies below 
1 keV, i.e. the data points in the spectra are shifted to 
the right. This causes the typical residual pattern seen in 
the bottom panel of Fig. 4 where data and model emis- 
sion lines are shifted slightly against each other in energy. 
This should affect the derived abundances only slightly 
and will not affect our conclusions. Other residuals near 
2.2 keV are caused by problems in the effective mirror 
areas around the Au edge. 

The modulation of the source is mainly due to a strong 
PL absorption varying from 13.0xl0 22 cm~ 2 in the dip to 
~ 6xl0 22 cm~ 2 in quiescence and down to (1.5 - 2.2) xlO 22 
cm~ 2 during flares. In the same time, the soft RS compo- 
nent is affected by a much lower (0.04 - 0.1) xlO 22 cm~ 2 
absorption increasing only to 0.2xl0 22 cm -2 during the 
dip. However, the flares are marked also by a significant 
increase in the RS intensity. 
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Table 3. Best- fitting spectral parameters 





Power-law 
Apl (*) a (*) 


N H 


Thermal component 
A RS kT (kcV) (*) N H 


XrVdof 


Flare 1 
Quiescent 
Flare 2 
Dip 


0.0230 (12) 1.35 (3) 


2.25 (11) 
6.43 (28) 
1.48 (10) 
13.1 (6) 


0.021 (3) 0.64 (5) 0.041 (6) 
0.009 (2) 0.145 (26) 
0.033 (4) 0.044 (6) 
0.011 (2) 0.22 (3) 


1.27/2467 



(*) single parameter in the joint fit to all four spectra 

Apl in photons cm -2 s" 1 keV" 1 at 1 keV, a photon index, Nh in units of 10 22 cm -2 
Ars emission measure in units of 10~ 14 /(4-7rd 2 ) J n e nHdV, 

where d is the distance to the source (cm), n e and nH are the electron and hydrogen densities (cm -3 ) 



5. Discussion 

The very good sensitivity of the EPIC cameras at energies 
below 2 keV allows for the first time to build a consistent 
picture of the source. A straightforward interpretation of 
its variability can be made with an unique effect of absorp- 
tion in front of a stable underlying emission with two sep- 
arated components. The hard component, best described 
by a power-law, is suggestive of a Comptonized-type spec- 
trum with a high energy break high above the 12 keV 
XMM-Newton limit. This same component was also iden- 
tified in EXOSAT and ASCA observations (PA86, Church 
et al. 1998). However at low energy, the XMM-Newton ob- 
servations show clearly that the emission is dominated by 
a thin hot plasma, affected by a much lower absorption 
than the hard component. This is at variance with previ- 
ous interpretations involving a second unabsorbed Comp- 
tonized spectrum (PA86) or a blackbody emission (Church 
et al. 1998). 

A constraint on the size of the emission area for 
the hard Comptonized emission is provided by its nearly 
complete eclipse in the orbital cycle. For a binary orbit 
with a P or b = 3.82 h, the transition time (ST) in and 
out of eclipses yields a typical size of R = 5.41 x 10 11 
(M x +M c )^ 3 -(cvT/P orb )-sim cm, where and M c are 
the masses of the compact star and of the companion in 
Mq and i the inclination of the system. Eclipse transition 
times are best determined by high time resolution data 
of 19 eclipses seen by EXOSAT (PA86) and 10 eclipses 
seen by RXTE (Hertz et al. 1997) which give mean val- 
ues of - 6 s (EXOSAT) and - 5 s (RXTE). Assuming a 
compact star of 1.0 Mq, a 0.45 Mq companion and an 
inclination of i = 75° (PA86), yields a typical radius of 
the Comptonized emission of R = (2.4±0.5)xl0 8 cm. 

On the other hand, the soft thermal component is 
seen un-eclipsed so that its size should be significantly 
greater than the companion occulting sphere which, for 
a Roche lobe filling star on the main sequence is of ra- 
dius R - 7.9xl0 9 -P orb (hr) cm (Frank et al. 1992). This 
corresponds to ~ 3xl0 10 cm for EXO 0748-676. 



The picture of the system as seen by XMM-Newton is 
therefore a rather compact (~ 2xl0 8 cm) Comptonized 
region as a central hard X-ray source surrounded by a 
more extended (~ 3x 10 10 cm) hot halo- type thermal emis- 
sion. This is the opposite to the previous descriptions in 
which the central source was expected to be a (~ 2 kcV) 
blackbody surrounded by a Comptonized emission from 
the accretion disc corona (ADC) (Church et al. 1998). We 
stress that the inner Comptonized region is compatible 
in size with the disturbed inner edge of the accretion disc 
where low frequency QPOs are expected to form. If the ob- 
served (0.58 - 2.44) Hz QPOs (Homan et al. 1999) arc in- 
deed produced at the Keplerian frequency in the accretion 
disc, their distance to the centre is given by r — 1.50 10 s 
i/ _2 / 3 -(M x /Mq) cm, where v is the QPO frequency, which 
corresponds to a range of (0.8 - 2.2) xlO 8 cm in radial dis- 
tance. The ADC is then strictly restricted to the central 
part of the disc. The spectral variability demonstrates that 
we always see this inner ADC strongly absorbed even dur- 
ing flares, most probably because of the high inclination 
and of the occultation by the external parts of the disc 
including the bulge during the dip. 

The most extended halo region is much less absorbed. 
In fact we note that the lowest best fitting value of Nh 
is a factor two lower than the average galactic column 
density in this direction Nh ~ l.lxlO 21 cm~ 2 (Dickey & 
Lockman 1998). This is possibly due to the more com- 
plex shape of the continuum in this region and a con- 
tamination by highly ionized regions in the system which 
should be accounted for by a warm absorber. The ab- 
sorption during the dip is a factor ~ 2 higher than the 
galactic value, revealing that the thermal extended halo is 
also significantly affected by absorption of the bulge. This 
is in favour of a flattened rather than spherical halo. The 
emission measure of this thermal halo in quiescent state 
is ~ 0.8 xlO 58 cm" 3 for a 10 kpc distance. For a spher- 
ical halo, with a typical dimension of ~ 3xl0 10 cm, the 
electron density will then ben~ 0.8xl0 13 cm" 3 . If a flat 
geometry is assumed, with a radius of ~ 3xl0 10 cm and a 
typical height of h ~ 10 9 cm, this density is significantly 
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Fig. 4. The XMM-Newton spectrum of EXO 0748-676. 
Spectra of four different parts of the orbit al light curve 
are shown such as flaring episodes Fl (red) and F2 (blue), 
inter- flare Q (green) and dip D (light blue). The spectra 
were multiplied by different factors to separate them in 
the plot (Fl 9.0, F2 27.0, Q 3.0, D 1.0). The best fit- 
ting model (histograms) is a combination of a highly ab- 
sorbed power-law and a thin plasma (0.6 keV) thermal 
emission component with non-solar abundances and sep- 
arate (smaller) absorption. The bottom panel shows the 
contribution of the residuals to the total \ 2 - The spectra 
are reproduced by varying mainly the amount of absorp- 
tion in the power-law component. Note that the intrinsic 
power-law intensity is constant, i.e. the spectra at the high 
energy end, when not scaled, lie on top of each other 



higher at n ~ 5.3 xlO 14 cm" 3 . We note that very simi- 
lar values of emission measures and densities are derived 
from the characteristics of the absorption edges and emis- 
sion lines of highly ionized elements clearly detected by the 
RGS at energies lower than 2.5 keV (Cottam et al. 2001). If 
the central source is seen through such a flat halo, a rough 
estimate of the absorption is Nh ~ h.n ~ 5.3x 10 22 cm" 2 . 
This value is consistent with the measured absorption of 
6.4xl0 22 cm" 2 (Table 3). The estimated dimension and 
density of the halo suggest that it may be a hot atmo- 
sphere at the surface of an irradiated accretion disc. The 
high inclination of the system will impose then that the 
central source is always seen through this extended highly 
ionized region on top of the accretion disc. A better de- 
scription of the halo may therefore be in terms of a fully 
photo-ionized model which is beyond the scope of this let- 



ter. Some caution should then be taken in interpreting the 
abundances derived from a pure collisional RS model. 

Assuming a distance of ~ 10 kpc, consistent with a 
maximum Eddington luminosity in the bursts (Gottwald 
et al. 1986), the (2 - 10 keV) unabsorbed luminosity of 
the source is constant at ~ 2.0xl0 36 erg s" 1 which is 
comparable to the luminosities observed by ASCA and 
EXOSAT, though at the lowest end. In the total XMM- 
Newton range, the (0.5 - 10 keV) luminosity only slightly 
varies from 2.64 (in quiescence) to 2.79xl0 36 erg s _1 
(flares). In this range, the thermal halo luminosity varies 
from 0.6 (quiescence) to 2.1 xlO 35 erg s" 1 (flares), con- 
tributing from 2 to 7.5% to the total luminosity. The in- 
crease of the halo luminosity is apparently not linked to 
any simultaneous changes in the central source luminosity. 
The total increase of energy inside the flares with typical 
duration of ~ 2000 s is ~ 3xl0 38 erg. This is only a small 
fraction of the energy ~ 2xl0 39 erg released during the 
bursts. It is then possible that the halo changing luminos- 
ity may be due to a delayed re-processing of a previous 
burst at the surface of the accretion disc. 
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